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Ultra-cold atoms which are subject to ultra-relativistic dynamics are investigated. By using 
optically induced gauge potentials we show that the dynamics of the atoms is governed by a Dirac 
type equation. To illustrate this we study the trembling motion of the centre of mass for an effective 
two level System, historically called Zitterbewegung. Its origin is described in detail, where in 
particular the role of the finite width of the atomic wave packets is seen to induce a damping of 
both the centre of mass dynamics and the dynamics of the populations of the two levels. 
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I. INTRODUCTION 

Ultra-cold atoms obeying ultra-relativistic dynamics 
would be a match made in heaven. An atomic cloud 
cooled down to nano-Kelvin temperatures offers an un- 
precedented opportunity to study and manipulate a true 
quantum gas. Relativistic dynamics, on the other hand, 
seems at first sight incompatible with the concept of an 
ultra-cold quantum gas. This is not necessarily the case. 
It has been noted in earlier works on atomic gases [H[2J[3], 
but also in layers of graphene [4] , that it is indeed possi- 
ble to study relativistic dynamics for Systems which are 
inherently non-relativistic. Counter-intuitively this is the 
case in the low-momentum limit. 

Already in the early days of quantum mechanics the 
dynamics of relativistic particles attracted a lot of at- 
tention. It was soon realised that a number of counter- 
intuitive results would follow in the relativistic limit, such 
as the Zitterbewegung, i.e., the trembling motion of the 
centre of mass of a wave packet. Zitterbewegung has at- 
tracted a lot of interest over the years and is continuing 
to be an active field of research which involves a broad 
ränge of physical Systems [H |6j |7] . 

The notion of Zitterbewegung itself has its roots in 
the work of Schrödinger on the motion of the free parti- 
cle based on Dirac 's relativistic generalization of a wave 
equation for spin-1/2 particles [8j|9j[lÜ]. From these early 
days on the existence of the Zitterbewegung for relativis- 
tic particles has also been the subject of some controversy 

hd na. 

Zitterbewegung results from the interefence of the pos- 
itive and negative energy Solutions of the free Dirac equa- 
tion. The frequency of this interference process is deter- 
mined by the energy gap between the two possible Solu- 
tion manifolds. In case of the historically first discussed 
motion of a free electron wave packet this energy gap is 
on the order of twice the rest energy 2m e c 2 of the elec- 
tron, i.e., the energy necessary to create an electron-hole 
pair. For a complete description the coneepts of quan- 
tum field theory, i.e., particle creation and annihilation, 



have to be introduced. If the Zitterbewegung is driven 
by a process which corresponds to an energy of the Or- 
der of the rest mass, then this unfortunately also makes 
it rather unlikely to observe the trembling motion with 
real electrons [T2"] . 

A number of physical Systems, can, however, be de- 
scribed by an effective Dirac equation. For these Systems 
the creation energy of all partieipating real particles is 
much larger than the gap energy and thus the processes 
of particle creation and annihilation can be disregarded. 
Moreover, the considered particles are not real but of a 
quasi-particle nature, e.g. spin-states. The combination 
of these features allow us to study the phenomenon of 
Zitterbewegung in regimes far away from its initial dis- 
covery - the motion of a free electron. 

In this paper we will study the trembling centre of 
mass motion of a two-level System which is subject to an 
off-diagonal matrix gauge potential. In the limit of low 
momenta and strong gauge fields the dynamics is well 
described by a Dirac- type equation. In recent papers [2, 
I13j it has been shown how atoms with an internal tripod 
level strueture, see Fig. [T] may evolve under the influence 
of an effective non-Abelian vector potential. Here we 
restrict the motion to only one direction (see also Ref. 
[3] for a two-dimensional description of Zitterbewegung) . 

The paper is organized as follows. In the following 
section we will briefly outline the derivation of the gauge 
potentials for the spin System. As an example of the re- 
sulting dynamics in the presence of non trivial gauge po- 
tentials we study the Zitterbewegung for neutral atoms. 
Finally we discuss the phenomenon of Rabi-type oscilla- 
tions occurring in this context and the damping mecha- 
nisms due to finite size effects. 



II. THE EQUATION OF MOTION 

In the following we will assume the motion of the atoms 
to be restricted to one dimension. We choose our coordi- 
nate System such that the x-axis is aligned along that par- 



2 



ticular dimension. A gas of ultracold atoms can be con- 
sidered dynamically one-dimensional if the corresponding 
transversal energy scale given by the transversal trapping 
frcqucncy is much higher than all other energy scales, 
such as the temperature or chemical potential in the pres- 
ence of collisional interactions. 

To the effectively one-dimensional cloud of cold atoms 
we apply the schcmc for inducing non-Abelian gauge Po- 
tentials as presented in [2], and obtain in the limit of 
low momenta a quasi-relativistic Situation as shown in 
[2]. For this purpose, we consider the adiabatic motion 
of atoms in the presence of three laser beams. The tech- 
nique is remarkably versatile and offers the possibility to 
shape the gauge potentials quite freely. Various possibil- 
ities exist for creating non-trivial equations of motions 
[13] . Here we have chosen a laser configuration wherc 
two of the beams have the same intensity but counter- 
propagate. The third laser beam has a different intensity 
compared to the two other laser beams. Its wave vector 
is chosen to be perpendicular to the axis defined by the 
propagation direction of laser 1 or 2. The configuration 
is depicted in Fig. [T] 

By defining the total Rabi-frequency £1 = 
and 



the mixing angle from 
| 2 / 1 3 1 , we can write the Rabi 



tanö = y |ni|- -i- [ü 2 
frequencies of the participating laser fields in the follow- 
ing form: Q x = ^ sin 0e~ iKX /V2, fl 2 = Q sin 0e tKX /V2 
and f^3 = cos 0e~~ lKy . Applying this notation we 
find that the interaction Hamiltonian after the dipolc- 
approximation and in the interaction picture is given 

by 



H; 



ft fii|o)(i| + n 2 |o)(2| + n 3 |o)(3| 



(i) 



in the case of resonant laser fields. The Hamiltonian 
Hi n t yields two dark states |-D,), i = 1,2, which contain 
no contribution from the excited State 10): 
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- iKV (e iKX \l) -e~ lKX \2)), 



(2) 
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Both dark states are eigenstates of li int with zero 
eigenenergy. They depend on position due to the spa- 
tial dependence of the Rabi frequencies f2j . 

The bright State \B) ~ + fl* 2 \2) + S7g|3> is cou- 

pled to the exited state |0) with the Rabi frequency Q 
and therefore separated from the dark states by energies 
±Ml. If |f2| is large enough compared to any two-photon 
detuning or Doppler shifts due to the atomic motion, we 
can neglect transitions out of the dark states, i.e., we use 
the adiabatic approximation. 

In this limit it is sufficient to expand the general state 



vector |x) of the quantum System in the dark state basis 

(4) 
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where ^(r) are the expansion coefficients. These coeffi- 
cients represent the centre of mass motion of the atoms 
in the dark state i. By collecting the wave functions in 
the spinor 



* = 



#2 



(5) 



we find that the latter obeys the effective Schrödinger 
equation 
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(6) 



where p x dcnotes the momentum along the x-axis and m 
is the atomic mass. Here A is an effective vector potential 
matrix, also called the Mead-Berry connection [T3 HS] 
and V and <& are effective scalar potentials matrices. The 
gauge potentials A„, m = ih(D n (r)\\7D m (r)) and $„, m = 
2^ J2i A n ,/A; iTO emerge due to the spatial dependence of 
the dark states. The other scalar potential is defined by 
V n . m = {D n {v)\V\D m {v)) with V = E"=i Vj{r)\j)(j\ and 
Vj(r) being the trapping potential for atoms in the bare 
state j. 

With the setup presented in Fig. [T] these potentials 
take in x-direction the form 
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(7) 
(8) 
(9) 



where we have introduced the notation n' = kcos(ö) and 
assumed that the external trapping potentials for the first 
two atomic states are the same, i.e. Vi = Vi. In addition, 
the external trap in the transversal direction is assumed 
to dominate over any effective gauge potential in this di- 
rection. Hence, we can use an effectively one-dimensional 
equation of motion. 



III. THE DIRAC LIMIT 

The existence of an energy gap is necessary in order to 
observe Zitterbewegung. Such a gap is obtained by the 
different, but now constant, trapping potentials V± and 
V3 which can be altered by detuning the corresponding 
lasers from the atomic transitions. Alternatively, the in- 
tensity ratio of the laser beams can be used to adjust the 
scalar potential. 
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FIG. 1: One possible laser configuration for the tripod Sys- 
tem which results in a non-trivial gauge potential for the two 
corresponding dark-states. 



It is convenient to introduce the following notation 
1 



V~ = 



Vu + $u - (Vi 



22 



22 



(10) 



and shift the zero level of energy. The trapping potential 
then reads 



V + $ = V z a z 



V z 



1 

-1 



(11) 



where a z is one of the Pauli spin matrices. In the limit 
of low momenta, i.e., \p\ <C Hk, we can neglect in Eq. 
([ö]) the kinetic energy term and are left with an effective 
Dirac equation 



[C<JxVi 



2m 



V z ha z 



(12) 
(13) 



Here 



Kk 



is an effective recoil velocity, which is typ- 



ically on the Order of cm/s for alkali atoms. The Dirac 
limit in Eq. (13) is most clearly justified by considering 



the corresponding linearised dispersion relation from Eq. 
^ [2]. Equation (13l is the starting point of our main 
discussion. We note that the A 2 can be absorbed into 
the potential term and we are hence left with an equation 
which resembles the Dirac equation for a free relativistic 
particle with the rest energy substituted by the potential 
energy difference between the two levels. 



IV. ZITTERBEWEGUNG 

In most textbooks [T7J [TS] Zitterbewegung is derived 
by solving the Heisenberg equation for the Space Oper- 
ator. For a free particle with rest mass m, the Dirac 
Hamiltonian Hu containing the speed of light c 



Ho = cap + ßmc 



(14) 



is used to obtain the time dependence for the Space Oper- 
ator x. By using the anticommutation properties of the 
Dirac a and ß matrices one finds 



*(*) 



x(0) 

ihc . 



cpffß 1 (15) 




FIG. 2: Density plot showing the Dirac limit and the role of 
the initial width of the wave packets. Figures (a)-(c) display 
the füll Schrödinger dynamics using eq.|6| for Gaussian initial 
states with increasing width er. This results in a sharper mo- 
mentum distribution with p <C hu increasingly fulfilled. The 
pure Dirac case is shown in figure (d) for comparison. The 
dynamics in (a)-(d) shows in addition to the Zitterbewegung 
also a damping. 



In equation ( 15 1 the first and second term describes a 



motion which is linear in time, while the third term gives 
an oscillating contribution, the Zitterbewegung. To ob- 
serve this trembling motion an initial 4-component spinor 
State needs to contain positive and negative parts as the 
a matrix is mixing these. The frequency of the oscillat- 
ing term can be estimated in the particles rest frame as 
2mc 2 /H. This typically large energy is the energy differ- 
ence between a particle and antiparticle. 

To obtain Zitterbewegung in cold quantum gases we 
emphasise that the two dark states described by equation 
( 13 1 are not a particlc-antiparticle pair, but they are still 



separated by an energy gap which is generated by the 
constant potential term in Eq. (11). Different to previous 



work [5] the Dirac Hamiltonian now contains an effective 
rest mass V z . 



V. DARK STATE DYNAMICS 

Our System with two degenerate dark states shows not 
only a relativistic behaviour, but also properties familiär 
from two level Systems in quantum optics |19j . 

In order to see this in more detail we write the spinor 
ty(x) as a combination of slowly varying envelopes, 4>i(x) 
i = 1,2, and coefficients which describe the population 
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FIG. 3: Left column: The density as a function of time shows 

fi 2 2 

Zitterbewegung for different energy gaps ((a) with V z = 

and (c) with V z = 3 fi 2l ^ ) and in (e) with an initial momen- 
tum fco = k' . Right column: The centre of mass shows the 
expected oscillation ((b) and (d)) with an upward drift. With 
a initial momentum the behaviour is different, as can be seen 
in (e) and (f), where the Zitterbewegung breaks down after 
a few oscillations as the two states are moving in different 
directions. The initial spinor was in (a)-(d) (1, l,) T /\/2 and 
in (e) and (f) (1, e" r/4 ) T /\ / 2- 




FIG. 4: Rabi type oscillations superimposing the Zitterbe- 
wegung can be seen for an initial preparation with nonzero 
momentum fco = K 1 . The dynamics is sensitive to the ini- 
tial conditions. Figure (a) and (b) show |^i(£,a;)| 2 for the 
initial spinor spinor (1, l) T /\/2 and (1, e I ' r ' /4 ) T /v / 2 respec- 
tively, whereas figure (c) and (d) show | \I/2 (i, x)\ 2 for the initial 
spinor (1, 1) T /V2 and (1, e i7r/4 ) T /V2. The population of the 
dark states is depicted in (e) with (1, l) T /\/2 and in (f) with 
(l,e ln/4 ) T /V2. 



of the two dark states, 
# = 



<t>i(x)ci(t) 

(ß2(x)c 2 (t) 



(16) 



The spatial shape <f>i(x) should change much slower in 
time than the population Ci(t) of the z:th component of 
the spinor. The Solutions are normalised according to 
{4>i\4>i) = 1 an( i l c i| 2 4- |e 2 | 2 = 1. After inserting the 



ansatz from equation (16l into equation (13) we obtain 



a set coupled differential equations for the coefficients: 



where 



and 



^ = -z{4>2\Px\(pl) 



(17) 



(18) 



(19) 



The two spin components are coupled, hence the Solu- 
tions to equation (17l will show oscillations between the 



two dark states with a frequency 



4 



im 



v z . 



(20) 



For a vanishing overlap integral f2, the coupling between 
the spin components in Eq. (17) becomes zero and we 



expect no oscillations of the populations. 

If the initial State has a nonzero momentum of the form 



cj>i{x) 



-x ja -\-ikox 



(21) 



with a width a and momentum fco, we obtain a non-zero 
f2 which will consequently result in population transfer 
between the two dark states. 
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VI. EXACT SOLUTIONS IN THE 
SCHRÖDINGER LIMIT 

The Zitterbewegung is typically a transient phe- 
nomenon in Systems such as graphene and carbon nan- 
otubes |2nj- The wave packets corresponding to the pos- 
itive and negative energies move in different directions 
which will inevitably cause any interference effect to de- 
cay. In the atomic case this is not necessarily so. With 
the coupled atom-light sctup thcrc is no attcnuation of 
the Zitterbewegung for wave packets in the Dirac limit 
without initial momentum, see Figs. [3](a)-(d). The pre- 
dicted attenuation [5Ü] is, however, observed if the initial 
states have a nonzero momentum, as shown in Figs. |3je) 
and[3jf). In this case, an attenuation of the amplitude of 
the dark State population diffcrence also occures. This is 
illustrated in Fig. |4j 

The Situation is typically more complicated due to the 
finite size of the atomic wave packet. The atomic cen- 
tre of mass motion is strictly speaking always governed 
by the Schrödinger equation ([6]), and only in the limit of 
small momenta compared to the corresponding momen- 
tum from the effective gauge potential, are we allowed to 
construct a Dirac type equation. This, however, would 
also require us to only use wave packets with zero mo- 
mentum spread. In the following we will analyse the role 
of the finite width of the wave packets. 

The exact Solution of the Schrödinger equation ^ with 
the gauge potentials from Equation ([7])-(|9|, can readily 
be written down in momentum space, 



*(fc,r) = e 



-i{k +2a m k+V z <T, 



> T *(M), (22) 



where now the dimensionless fc is expressed in units of 
k and the time r in units of 2m /hn 2 . If we choose a 
Gaussian momentum distribution with a width A for the 
initial State, 



To illustrate this we choose an initial State with c\ — 
C2 = \j\[2. The resulting populations of the two dark 
states are then given by 



hMr)| 2 = 

l^ 2 (r)| 2 : 



Ä|*i(0,/c)| 2 (l 

i 

dfc|* 2 (0,fc)| 2 (l 



2kV z 



2Wy 



2 Kr)) (26) 
\üü kT )){27) 



The population difference, AN(t) = \ipi(t)\ 2 - IM*)) 2 , 
can easily be calculated in the limit A = 0. In this limit 
the Gaussian initial State turns into a representation of 
the delta funetion. 



AiV(r) 



lim 

A^O 



dk- — =e 



sin 2 (w fc r) 



4k Q V z 



4k 2 



:sin 2 ( 



V? 



4k 2 + V 2 r)m 



From this rcsult we see the importance of the initial mo- 
mentum fco- For fco = there is no transfer of popu- 
lation between the dark states, whereas for a nonzero 
initial momentum the amplitude of the population oscil- 
lation is proportional to fco. In addition, the frequency 

4fcg + V 2 is fco dependent as well. 
The expectation value for the centre of mass can be 
calculated from 



(x(t)) 



dfc* t (fc,r)9 fc *(fc,r 
4fc 2 



dke a 2 



V 2 

^-sin(2w fc r))(29) 



From the first term under the integral sign we obtain a 
drift term for the centre of mass, 



*(fc,0) = -1 e -C*-*o) 2 / 2 A* ( Cl \ ( *i(M) 

(23) 

we obtain an exact time dependent Solution of the form 



ci cos(w fe r) + ' (ciy "t C22fc) sin ( w fc T ) 



, (24) 



\ ( \ iiczV z — ci 2k) • / \ 

\ c 2 cos{uJkT) - J '— >-sm(u; k T) 

where we have introduced the fc-dependent frequency 



4fc 2 + V 2 



(25) 



With the Solution of ^>(k,r) we can calculate the centre 
of mass motion for the two-component wave packet where 
we use the Standard definition of the density, p(k,r) — 
|\E'i(fc,r)| 2 + | v E , 2(fc J T )\ 2 1 or indeed any dynamics in the 
populations of the two dark states. 



Xd 



1 



' tt -4- e Erf ^c ( — ) 



A 



A 



(30) 



where Erfc is the complcmcntary Error funetion. In the 
limit of V z /A 3> 1, and using the asymptotic expansion 
of the Error funetion, 



Erfc(x) 



n 



n\(2xY 



(31) 



we obtain a reduced drift as a funetion of increasing 
V z /A. This is a finite size effect and stems from the 
finite width of the wave packet. 



The second term under the integral in ( 29 1 Stands for 



the Zitterbewegung. Again, in the limit V z /A 3> 1, the 
integral can be readily calculated and gives for the oscil- 
lating part 



^ sin \ 2V z t + | arctan(^r) 



V z 



(1 



A 4 
Kil'- 



2)1/4 



(32) 
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Prom this expression we see that a spread in the mo- 
mentum distribution will cause a damping also for the 
oscillating term of the centre of mass. The damping of 
the Zitterbewegung is relatively slow, but is inevitable. 
The underlying equation is after all the Schrödinger equa- 
tion, and only in the limit A — can we strictly spcaking 
use a Dirac-type equation which results in no expansion 
of the wave packet, hence no damping. With the füll 
Schrödinger equation a free wave packet will always ex- 
pand, albeit slowly if A is small, and hence will also show 
a damped Zitterbewegung. 

For a typical alkali atom such as 87 Rb with a wave 
packet width of 10 /im onc would get A 2 t/AV z > 1 for 
times larger than 1ms, with a centre of mass oscillation 
frequency of the Order of IkHz. Hence a broad wave 
packet as initial state would favour the detection of the 
Zitterbewegung. 

VII. CONCLUSIONS 

In this paper we have showed using ultra-cold atoms 
how Zitterbewegung, known from relativistic physics, is a 
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generic phenomena which will naturally occur in Systems 
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in a position to for instance study a System which would 
correspond to a confined Dirac particle by introducing ex- 
ternal atomic potentials, either by optical or magnetical 
means, for the atoms. This consequently leads us to pon- 
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atoms will play an important role, where one would be 
faced with a nonlinear Dirac-type equation to describe 
the dynamics. The source of the nonlinearity is however 
nontrivial due to the underlying collisions between the 
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